Dimensionality selection in a molecule-based magnet 
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Gaining control of the building blocks of magnetic materials and thereby achieving particular 
characteristics will make possible the design and growth of bespoke magnetic devices. While progress 
in the synthesis of molecular materials, and especially coordination polymers, represents a significant 
step towards this goal, the ability to tune the magnetic interactions within a particular framework 
remains in its infancy. Here we demonstrate a chemical method which achieves dimensionality 
selection via preferential inhibition of the magnetic exchange in an S = 1/2 antiferromagnet along 
one crystal direction, switching the system from being quasi-two- to quasi-one-dimensional while 
effectively maintaining the nearest-neighbour coupling strength. 



Coordination polymers are self-organising materials 
consisting of arrays of metal ions linked via molecu- 
lar ligands, with non-coordinated counterions supplying 
charge neutrality. The choice of initial components per- 
mits a high level of control over the final product, en- 
abling many different polymeric architectures to be ob- 
tained pQ . These materials provide a route to successful 
crystal engineering, and a number of functionalities are 
being actively studied, including gas storage [2HI]> opto- 
electronic [5] [B] , ferroelectric [7J [5] and magnetic proper- 
ties [hhh]- 

Although it is now possible to generate an assortment 
of disparate magnetic lattices using this method [TSJ US] , 
true control of magnetic exchange interactions implies an 
ability to adjust selected parameters while keeping oth- 
ers constant. To this end, a series of coordination poly- 
mers based on Cu(II) ions bridged by pyrazine (C4H4N2) 
molecules have proven to be highly versatile. In these 
systems it has been shown that it is possible to alter 
significantly the primary exchange energies via adjust- 
ment of the ligands [T7J and the counterions [TS1 IS9j , or 
fine-tune the exchange by a few percent via isotopic sub- 
stitution [20] . all the while maintaining the same basic 
metal-pyrazine network. In this paper we demonstrate 
the power of this strategy by chemically engineering a re- 
duction in the dimensionality of a magnetic system. After 
first designing a material based on coordinated planes of 
Cu(II), we adapt the recipe such that the ligand bridges 
are broken along a specific crystal direction, resulting 
in a chain-like compound. Because the ligand mediating 
the magnetic interactions in both cases is unchanged, the 
nearest-neighbour exchange energies of the two materi- 
als are found to be equal to each other to within 5%. 
The difference in numbers of nearest-neighbours, how- 
ever, means that the strength of the combined exchange 



interactions acting on each magnetic ion in the quasi- 
two-dimensional material is twice that of its quasi-one- 
dimensional cousin. 

Figs, [ija) and (b) show the crystal structure of 
orthorhombic [Cu(pyz)2(pyO)2](PFg)2 (where pyz = 
pyrazine and pyO = pyridine-iV-oxide, C5H5NO) deter- 
mined using single-crystal x-ray diffraction [21] . S = 1/2 
Cu ions are linked by pyz molecules into nearly square 
planar arrays, with perpendicular non-bridging pyO lig- 
ands keeping the planes well-separated. Because of the 
separation, as well as the staggered arrangement of adja- 
cent planes shown in Fig.[IJb), magnetic exchange ener- 
gies are likely to be very small along the c-direction. In 
contrast, Cu — pyz — Cu bridges are known to be good me- 
diators of antiferromagnetic superexchange [221 123j and 
so the magnetic properties of this material are expected 
to be quasi-two-dimensional. This is confirmed by the 
magnetic measurements described below. Sample syn- 
thesis involves mixing together of the molecular compo- 
nents in a solution of water and ethanol. Intermolecular 
self-organisation means that only a small amount of inter- 
vention is subsequently required. To achieve the desired 
planar structure the pyz and pyO molecules were added 
in a 3:1 ratio, previous experience suggesting that to ac- 
count for the potential for pyO to substitute for pyz, the 
ligands must initially be in a proportion different to that 
found in the final product. In order to create a similar 
sample, but one based on Cu — pyz chains rather than 
planes, we reduce the pyz:pyO ratio to 2:1 and proceed 
with the synthesis in a similar way. The resulting mate- 
rial has the composition [Cu(pyz)(pyO)2(H 2 0)2](PF 6 )2 
and the structure is shown in Figs.JIJc) and (d). Here the 
pyz ligands link Cu ions along the b-axis only, the other 
ligands being non-bridging pyO and water molecules. 
The alteration in composition has the effect of chang- 
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FIG. 1. (a) View of the crystal structure of the planar material [Cu(pyz)2(pyO)2](PFe)2 determined using x-ray diffraction 
showing the 2D Cu — pyrazine network in the afo-plane. (b) a projection of the same structure along the a-axis highlighting the 
shift between adjacent Cu — pyrazine layers and the arrangement of the PFs counterions. (c) Crystal structure of the chain-like 
material [Cu(pyz)(pyO)2(H20)2](PFe)2 showing the ID Cu — pyrazine chains in the afc-plane, and (d) a projection along the 
chains showing the arrangement of the pyO and H2O ligands and the PFe counterions. The dashed lines in (c) and (d) enclose 
one unit cell. Cu = brown, C = grey, N = blue, O = red, P = orange, F = green. Hydrogens other than those in the water 
molecules have been omitted for clarity. 



ing the symmetry of the crystal from orthorhombic to 
monoclinic and reducing the number of formula units in 
the unit cell, but, most importantly for the magnetic ex- 
change, the Cu — pyz — Cu linkages along the a-axis are 
removed without altering the b-axis Cu-Cu separation 
by more than a fraction of a percent (6.914 ±0.001 A for 
the planar material and 6.851 ± 0.001 A for the chain 
compound) . 

Recent heat capacity measurements on the planar ma- 
terial in zero magnetic field see no evidence of a mag- 
netic transition down to the lowest temperatures mea- 
sured [24] . However, thermodynamic probes are known 
to be less sensitive to transitions driven by interplanar 
couplings Sll while local probes such as muon-spin re- 
laxation (/i + SR) are much more effective at determining 
the antiferromagnetic transition temperature, Tm 122] . 
Our /i + SR data on this compound, shown in Fig. p2[a), 
exhibits a clear precession signal which develops below 
Tn = 1.71 ± 0.02 K [56], demonstrating long-range mag- 
netic order throughout the bulk of the sample. In con- 
trast, the /i + SR data for the chain-like compound exhibit 
no resolvable oscillations, see Fig.[2]jb), probably due to a 
smaller ordered moment, but can be fitted to the expres- 
sion A{t) = Aoe~ xt + An , where An represents the non- 
relaxing part of the signal. Both An and A rise markedly 
below T = 0.26 K due to crossover from a regime in 
which the relaxation is dominated by dynamic magnetic 
fluctuations to one dominated by quasistatic magnetic 
order, see Fig. j2^b) insets. From these fits we estimate 
T N = 0.27 ± 0.01 K for this material [2"T] . 

The type of magnetic order displayed by the two com- 
pounds can be deduced from their low-temperature mag- 
netization (see Fig. [3ja)). The form of our pulsed-field 



magnetization data up to saturation is in keeping with 
that expected for low-dimensional antiferromagnets. The 
slightly concave curve exhibited by the planar material 
is typical of quasi-two-dimensional (Q2D) antiferromag- 
netic interactions [18) , while the more extreme curvature 
shown by data from the chain-like sample is indicative of 
a quasi-one-dimensional (Q1D) magnetic lattice, where, 
for an ideal system, dM/dB is known to diverge at the 
saturation field [55]. To support these observations we 
compare the data with the results of low-temperature 
quantum Monte Carlo (QMC) simulations based on the 
Hamiltonian 

h = j ^2 s i ■ s 3 - + j± 12 s * • s j - 9» bB 12 S i- w 



(i ,j) II 



Here, for a Q1D (or Q2D) system, J is the strength of the 
exchange coupling within the magnetic chains (planes), 
Jl is the coupling between chains (planes), and the first 
and second summations refer to summing over unique 
pairs of nearest neighbours parallel and perpendicular 
to the chain (plane), respectively. Comparisons of the- 
ory and data are shown in Fig. 3(b). Theoretical mag- 
netisation curves were calculated in finite steps of the 
J±/J parameter for both Q1D and Q2D magnetic lat- 
tices [HI [21] . For both materials curves corresponding 
to the quoted values of J±/J gave the best matches with 
experimental data; the Q2D curve with J±/J = 0.00 for 
the planar sample, and the Q1D curve with J±/J = 0.02 
for the chain-like material. The predicted curve for a 
three-dimensional antiferromagnet is also shown for con- 
trast. The deviation of the data from the simulations 
close to the saturation field is likely due to the finite 
temperatures at which the experiments were performed. 
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FIG. 2. (a) Example muon-spin relaxation (/i + SR) spectra 
measured on planar [Cu(pyz)2(pyO)2](PFe)2. Data at dif- 
ferent temperatures are offset for clarity. Inset: Evolution 
of the precession frequency v with temperature. Long-range 
magnetic order is observed below 1.71 K. (b) /i + SR spectra 
measured on chain-like [Cu(pyz)(pyO)2(H20)2](PFe)2. In- 
sets: The evolution of the amplitude An (left) and relaxation 
rate A (right) with temperature indicates the onset of long- 
range magnetic order below about 0.27 K. Red lines are fits 
to functions described in the text and 1211. 



The saturation field (B c ) can be extracted from the 
pulsed-field data and is found to be 23.7 ± 0.8 T for 
the planar material and 12.8 ± 0.4 T for the chain-like 
material. At saturation, the Hamiltonian above implies 
that g/iB-B c = nJ + n±J±, where, for each spin, n is 
the number of exchange bonds of interaction strength 
J (n = 4 for Q2D, 2 for Q1D), and n± is the number 
of exchange bonds of interaction strength J± (nj_ = 2 
for Q2D, 4 for Q1D) [HI HOI El] . From this, by assum- 
ing that J±/J <C 1 and using the appropriate values 
for the ^-factor determined from electron-spin resonance, 
we estimate the primary exchange couplings, J, to be 
8.1 ± 0.3 K for the planar compound and 8.8 ± 0.2 K for 
the chain-like compound. 

Another estimate of the magnetic dimensionality 
comes from the temperature dependence of the low-field 
magnetic susceptibility. Fits of such data for both mate- 
rials results in estimates of J that are in accord with those 
derived above from pulsed-field magnetization data [21] . 
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FIG. 3. (a) Normalized pulsed-field magnetization of pla- 
nar [Cu(pyz)2(pyO)2](PF6)2 at T = 1.5 K, and chain-like 
[Cu(pyz)(pyO)2(H 2 0)2](PF 6 )2 at T = 0.5 K. (b) The results 
of quantum Monte Carlo (QMC) simulations of the low tem- 
perature magnetization for 3D (J±/J = 1), 2D (J±/J = 0), 
and ID (J±/J ~ 0.02) antiferromagnets. The lines are the 
pulsed-field data scaled by the saturation field (-Be), (c) 
The relation between the exchange anisotropy and the ra- 
tio of critical temperature and primary exchange energy in 
Q1D and Q2D antiferromagnets deduced from QMC simula- 
tions [28]. The circles indicate the materials reported here, 
(d) Anisotropy of the gr-factor in the planar and (e) chain-like 
compounds measured using ESR at 10 K and 1.5K, respec- 
tively. Red lines are fits to g{6) — (g^, y sin 2 8 + g\ cos 2 #) 5 . 



The relative sizes of Tn and J are indicative of 
the anisotropy of the exchange interactions in a low- 
dimensional magnetic system. Using quantum Monte 
Carlo (QMC) calculations Yasuda et al. [15] developed 
empirical relations between these values for S = 1/2 
Heisenberg antiferromagnets, the results of which are 
shown in Fig. [3jc) . For compounds considered here the 
exchange anisotropies are found to be \J±/J\ ~ 10~ 4 in 
the planar material and 10~ 2 in the chain-like material. 
These values are in keeping with the comparison between 
data and simulation shown in Fig. [3jb) . 

The magnetic lattice in a low-dimensional system is not 
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always obvious from an inspection of the crystal struc- 
ture [17l [29] . Verification of the equivalence of the mag- 
netic and structural planes and chains in our materials 
comes from the angle-dependence of the g-factor as de- 
termined by electron-spin resonance. To first approxi- 
mation, the Cu(II) ions in the planar material have a 
local octahedral symmetry with a tetragonal distortion 
along the c-axis. In such situations, the unpaired spin 
in a d 9 configuration is expected to occupy the d x 2_ y 2 
orbital [3D], and the g-factor parallel to the distorted z- 
axis takes a larger value than those in the a;y-plane (see 
discussion in Ref. 1^9")) . Accordingly, the data in Fig. ^d) 
show that the g-factor is lowest when the excitation field 
is applied in the a6-plane, implying that the direction of 
highest electronic orbital overlap is the Cu — pyz direc- 
tions, with the strong Q2D exchange interactions being 
mediated via the molecular orbitals of the pyrazine. For 
the chain-like material, the distortion is more compli- 
cated due to the lower crystal symmetry. Nevertheless, 
using two rotations, the smallest values of the g-factor 
are found to be along the Cu — pyz and Cu — pyO bonds, 
as shown in Fig. [3^e) , implying that again the magnetic 
orbital is d x 2_ y 2. This allows for the possibility of good 
exchange coupling along those bonds, but as the pyO 
molecules are non-bridging ligands, the Q1D magnetism 
must be mediated along the Cu — pyz chains. 

Taken together, these experimental observations paint 
a complete picture of the two closely-related magnetic 
systems. From the point of view of magnetic superex- 
change, the chain direction in the Q1D material looks 
very similar to the two Cu — pyz directions in the Q2D 
material and so, despite the compositional differences, 
the primary nearest-neighbour exchange energies remain 
largely unaltered. At the same time the critical field 
in the Q2D material is approximately double that for 
the Q1D compound because it has twice the number 
of nearest neighbours. The Q2D compound is strongly 
anisotropic, comparable to the most 2D materials yet 
identified [S6]. As suggested above, this is likely due to 
the disconnect and staggering that occurs between suc- 
cessive planes. The extreme anisotropy explains why the 
zero-field heat capacity was not sensitive to the antifer- 
romagnetic transition observed using /x+SR. In the Q1D 
material the anisotropy is less pronounced (even though 
the Tn/J ratio is smaller) because the chains are not 
staggered, there are twice as many next-nearest neigh- 
bours, and the shortest distance between chains (along 
the a-axis) is approximately half the interlayer separation 
in the Q2D material. 

These materials showcase the ability to take deliberate 
control over the magnetic properties of polymeric sys- 
tems. The self-organisation of the coordination polymers 
enables them to spontaneously form crystalline lattices 
whose structure can be anticipated with a high level of 
predictability. It is this predictive power, together with 
the ability to choose the starting ingredients and the 



knowledge of the exchange efficiency of various ligands 
accrued over the past few decades, that permits the pre- 
selection of exchange anisotropy. In the Cu — pyz systems 
we have demonstrated the ability to make this preselec- 
tion without significantly perturbing the magnitude of 
the primary interaction strengths, while previous stud- 
ies of similar materials have highlighted the capacity to 
tune the exchange couplings without changing the over- 
all dimensionality. Thus these compounds represent a 
promising approach to magnetic crystal engineering, and 
in particular raise the possibility of generating systems 
that exhibit higher ordering temperatures and other co- 
operative phenomena. 
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SUPPLEMENTAL MATERIAL ACCOMPANYING 
DIMENSIONALITY SELECTION IN A 
MOLECULE-BASED MAGNET 

Synthesis of planar [Cu(pyz) 2 (pyO) 2 ](PF 6 ) 2 . In 

the minimum amount of water, CuCl 2 -2H 2 (0.502 g, 
2.95 mmol) and AgPFg (1.489 g, 5.89 mmol) were mixed 
together to yield a white precipitate of AgCl and a pale 
blue solution containing the "Cu(PF 6 ) 2 " product. Tore- 
cover the blue solution, the AgCl solid was removed by 
careful vacuum filtration. In a separate solution, pyrazine 
and pyridine- A^-oxide in a 3:1 ratio were dissolved in a 1:1 
mixture of H 2 and ethanol to give a colourless solution. 
To this latter solution was added the "Cu(PF 6 ) 2 " solu- 
tion. Upon slow mixing of the chemical reagents, a small 
amount of pale blue powder formed, but was removed by 
vacuum filtration to afford a blue solution. Slow evapo- 
ration of this solution yielded a large mass of deep green 
square plates. 

Synthesis of the chain-like compound 
[Cu(pyz)(pyO) 2 (H 2 0) 2 ](PF 6 ) 2 . Using a reac- 
tion scheme similar to that above, the ratio of of 
pyrazine to pyridine- Af-oxide was decreased from 3:1 
to 2:1. Mixing of the chemical reagents immediately 
afforded a bright turquoise blue solution. The solution 
was covered with a piece of perforated Al-foil to allow 
slow solvent evaporation. Upon standing at room 
temperature for several weeks, small green blocks were 
deposited on the bottom of the beaker and subsequently 
isolated via vacuum filtration. 

Structural determination. Crystals of each com- 
pound were attached to a glass fibre and data were 
collected at 297(2) K using a Bruker/Siemens SMART 
APEX instrument (Mo K Q radiation, A = 0.71073A). 
Data were measured using omega scans of 0.3° per 
frame for 10 seconds, and a full sphere of data 
was collected. Data were analysed and cell param- 
eters retrieved using the supplied software. The 
data for [Cu(pyz)(pyO) 2 (H 2 0) 2 ](PF 6 ) 2 were found to 
be rotationally twinned and were deconvolved using 
CELLJMOW [Si] software to give a refined twinning ratio 
of 0.04±0.01. Absorption corrections were applied using 
TWINABS S2] . No such twinning was observed from the 
data for [Cu(pyz) 2 (pyO) 2 ](PF 6 ) 2 . Each cell component 
was refined using SAINTPlus [S3] for all observed reflec- 
tions. Data reduction and correction for Lp and decay 
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TABLE SI. Room temperature lattice parameters deter- 
mined from x-ray crystallography of [Cu(pyz)2(py0)2](PFe)2 
(Q2D) and [Cu(pyz)(pyO)2(H 2 0)2](PF 6 ) 2 (Q1D). *Note that 
this value is the perpendicular distance between Cu ions. 
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Cul-N8 (A) 
Cul-Nll (A) 
Cul-N14 (A) 

Cui-oi (A) 

Cul-02 (A) 


2.0452 ± 0.0019 
2.0646 ± 0.0024 
2.0690 ± 0.0024 
2.3164 ±0.0019 


2.044 ± 0.004 

1.969 ± 0.003 
2.420 ± 0.005 



TABLE S2. Coordination sphere bond lengths of 
[Cu(pyz)2(py0)2](PF 6 ) 2 and [Cu(pyz)(py0)2(H20)2](PF 6 )2 
(Q2D and Q1D, respectively) determined from room temper- 
ature x-ray crystallography. Atom numbering schemes can be 
found in Supplementary Figures SI and S2. 



were performed. The structure was found by direct meth- 
ods and refined by a least squares method on F2 using 
the SHELXTL program suite |S4j . The structures were 
solved in space groups P2\/m (#11) and Cmca (#64) by 
analysis of systematic absences. All non-hydrogen atoms 
were refined anisotropically. No decomposition was ob- 
served during data collection. 

The lattice parameters so determined are shown in Ta- 
ble IS 1 1 Also shown are the distances between Cu ions 
along the crystal axes. Alternate tilting of the pyz and 
pyO rings along certain directions, as well as stagger- 
ing of Cu — pyz planes (see Figs. 1 and 2), leads to unit 
cells that contain more than one formula unit. Table [S2l 
contains the coordination sphere bond lengths for both 
materials. The atom numbering schemes can be found in 
Supplemental Figures SI and S2, which respectively show 
the thermal ellipsoid plots for the planar and chain-like 
materials. 

Crystallographic data (.cif format) for the materials 
described here will be uploaded to the Cambridge Crys- 
tallographic Data Centre. 

Muon-spin relaxation. Zero field muon-spin relax- 
ation (/i + SR) measurements were made using the LTF 
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FIG. SI. Atom labelling scheme for the asymmetric unit 
of planar [Cu(pyz)2(pyO)2](PFe)2. Thermal ellipsoids are 
drawn at the 50% probability level. Symmetry transforma- 
tions used to generate equivalent atoms are #: -x, y, z; $: -x, 
y+1, z; %: -x+1/2, y, -z+1/2; &: x, y+1/2, -z+1/2. 




FIG. S2. Atom labelling scheme for the asymmetric unit of 
[Cu(pyz)(pyO)2(H20)2](PFe)2. Thermal ellipsoids are drawn 
at the 50% probability level. Only one orientation of the dis- 
ordered PF6 counterion associated with P2 is shown. Sym- 
metry transformations used to generate equivalent atoms are 
#: x, -y+1/2, z; $: x, -y+3/2, z; *: -x+1, -y+1, -z+1. 



instrument at the Swiss Muon Source at the Paul Scher- 
rer Institute. Powder samples were mounted in vacuum 
grease on a silver plate, which was fixed to the cold-finger 
of a dilution refrigerator. 

In a typical /i + SR experiment |S5| spin-polarized pos- 
itive muons are stopped in a target sample, where the 
muon usually occupies an interstitial position in the crys- 
tal. The observed property in the experiment is the time- 
evolution of the muon-spin polarization, the behaviour of 
which depends on the local magnetic field at the muon 
site. Each muon decays, with an average lifetime of 
2.2 fis, into two neutrinos and a positron, the latter par- 



ticle being emitted preferentially along the instantaneous 
direction of the muon spin. Recording the time depen- 
dence of the positron emission directions therefore allows 
the determination of the spin-polarization of the ensem- 
ble of muons. In our experiments positrons are detected 
by detectors placed forward (F) and backward (B) of the 
initial muon polarization direction. Histograms N-p(t) 
and JVb (t) record the number of positrons detected in the 
two detectors as a function of time following the muon im- 
plantation. The quantity of interest is the decay positron 
asymmetry function, defined as 



A(t) 



N F (t) 



N F (t) 



a, 



cxp 



N B (t) ' 



(SI) 



where a exp is an experimental calibration constant. A(t) 
is proportional to the spin polarization of the muon en- 
semble, and as a result depends upon the local magnetic 
field at the muon sites. 

As described in the main text, the planar sam- 
ple [Cu(pyz) 2 (pyO)2](PF 6 ) 2 exhibits oscillations in the 
/i + SR asymmetry below about 1.7 K, strongly indica- 
tive of long range magnetic order throughout the bulk 
of the material. This is because the local field causes a 
coherent precession of the spins of those muons for which 
a component of their spin polarization lies perpendicu- 
lar to this local field (expected to be 2/3 of the total 
spin polarization for a powder sample). The frequency 
of the oscillations is given by Vi = "f^Bi/2'Tr, where 7^ is 
the muon gyromagnetic ratio (= 27r x 135.5 MHz T _1 ) 
and Bi is the average magnitude of the local magnetic 
field at the zth muon site. The precession frequencies are 
proportional to the magnetic order parameter. The evo- 
lution of the muon precession frequency is shown in the 
inset to Fig. 2A of the main text and fitting this to the 
phenomenological function v{T) = u(Q) [1 - (T/T N ) Q f 
yields a transition temperature of Tn = 1.71 ± 0.02 K 
and exponents a = 1.1 ± 0.3 and = 0.22 ± 0.02. The 
analysis of /i + SR data for planar coordination polymers 
is described fully in Ref [S6] , 

The /i + SR data for the chain-like compound 
[Cu(pyz)(pyO)2(H 2 0)2](PF 6 )2 exhibit no resolvable os- 
cillations, but a monotonic relaxation of the muon polar- 
ization across the measured temperature regime (main 
text, Fig.2B). As previously mentioned, the data were 
fitted to the expression A(t) = A e~ xt + An, where A\\ 
represents the non-relaxing part of the signal. Around a 
temperature of T = 0.26 K we observe a large change in 
the behaviour of both Ay and A which we attribute to the 
onset of long-range magnetic order below that tempera- 
ture. The increase seen in Au on cooling is character- 
istic of a transition from a region of dynamic magnetic 
fluctuations to one of quasi-static magnetic order with 
decreasing temperature [S7]. This is attributable to the 
fact that above Tn dynamic field fluctuations will relax 
all muon spins; but upon ordering in a polycrystalline 
material, approximately 1/3 of the muon spins will be 
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aligned along the direction of the quasi-static local field 
and therefore will not be relaxed. These spins make a 
non-relaxing contribution to the spectra, detected as an 
increase in A\\ [S7 . In addition, the relaxation rate A 
is expected to vary as a simple function of the second 
moment of the magnetic field distribution (B 2 ) and its 
sudden increase below Tn reflects the increase in the size 
of the ordered moment (and hence (B 2 )) with decreasing 
T, below the transition. The onset of both effects occurs 
at the same temperature, from which we estimate the 
transition temperature of Tn = 0.27 ± 0.01 K. 

Hamiltonian and exchange constants. Our Hamil- 
tonian (main text, Equation 1) employs the single- J con- 
vention such that, in the case of just two S = 1/2 in- 
teracting spins, the separation between the singlet and 
triplet energies is equal to the interaction strength J. 
In some previous works, the Heisenberg model has been 
defined such that the singlet-triplet gap is 2 J [S8l [S9] . 
The exchange couplings derived from this definition are 
half those obtained by our model. Note also that, for 
materials with two different interchain exchange paths, 
J± in this Hamiltonian represents an effective interchain 
exchange. 

Pulsed-field magnetization. Measurements up to 
60 T (rise-time to full field ~ 10 ms) were performed 
at the National High Magnetic Field Laboratory in Los 
Alamos. Single crystals were mounted in 1.3 mm diam- 
eter PCTFE ampoules that can be moved into and out 
of a 1500-turn, 1.5 mm bore, 1.5 mm long compensated- 
coil susceptometer, constructed from 25 /xm high-purity 
copper wire. When the sample is within the coil and 
the field pulsed the voltage induced in the coil is propor- 
tional to the rate of change of magnetization with time. 
Accurate values of the magnetization are then obtained 
by subtraction of the signal from that taken using an 
empty coil under the same conditions, followed by nu- 
merical integration. The magnetic field is measured via 
the signal induced an adjacent empty 10-turn coil and 
calibrated via observation of dc Haas-van Alphen oscil- 
lations arising from the copper coils of the susceptometer. 
The susceptometer is placed inside a 3 He cryostat, which 
can attain temperatures as low as 0.5 K. Several orien- 
tations of single-crystal samples were used. The differ- 
ences between different orientations of the same sample 
were found to be small, and in keeping with the g-factor 
anisotropy determined by electron-spin resonance. The 
data shown in the paper were taken with B || a6-plane 
and T = 1.5 K for the planar material, and B || b and 
T = 0.5 K for the chain-like material. 

The magnitude of J extracted from the saturation field 
(taken to be the mid-point of the transition in dM/dB) 
explains why the pulsed-field data for the chain-like ma- 
terial is reminiscent of a Q1D antiferromagnet even at 
temperatures above Tn- The short-range correlations 
that occur above the transition give rise to a hump in 





Q2D 


Q1D 


technique 


Tn (K) 


1.71 ± 0.02 


0.27 ±0.01 


M+SR 


Be (T) 


23.7 ±0.8 


12.8 ±0.4 


M(B) 


J(K) 


8.1 ±0.3 


8.8 ±0.2 


M(B) 


J(K) 


8.10 ±0.01 


8.58 ±0.01 


X(T) 


9&y 


2.04 ± 0.01 


2.06 ±0.01 


EPR 


9z 


2.26 ± 0.01 


2.32 ±0.01 


EPR 


J±/J 


1 x 10~ 4 


1 x 10~ 2 


QMC 



TABLE S3. Magnetic properties of [Cu(pyz) 2 (pyO) 2 ](PF 6 )2 
(Q2D) and [Cu(pyz)(pyO)2(H 2 0)2](PF 6 ) 2 (Q1D) determined 
using ^i + SR - muon-spin rotation, M(B) - pulsed-field mag- 
netization, x(T) - low-field magnetic susceptibility, ESR - 
electron-spin resonance, and QMC - quantum Monte Carlo 
calculations. For the Q2D (Q1D) material the g-factors g xy 
and g z are measured with the magnetic field within the planes 
(chains), and parallel to the c-axis (perpendicular to the 
Cu — pyO bond), respectively. 



both the susceptibility and heat capacity Sll) of a Q1D 
Heisenberg antiferromagnet at temperatures T ~ J. The 
temperature at which the pulsed-field data were taken 
(T = 0.5 K) is considerably lower than this, and much 
closer to Tn = 0.27 K, implying that the correlations 
have grown to such an extent that for a non-local probe 
it is hard to distinguish this state from one with true 
long-range order. 

Quantum Monte Carlo simulations. The numerical 
results were obtained using the Stochastic Series Expan- 
sion (SSE) method to simulate the Hamiltonian of Equa- 
tion 1 on finite-sized lattices. SSE is a finite temperature 
QMC method based on the evaluation of the diagonal 
matrix elements of the density matrix in a suitable basis. 
The method is explained in detail in Ref. |Sllj and ref- 
erences therein. The simulations were performed on the 
NERSC computing facility in Berkeley, California. 

Temperature-dependent susceptibility. Measure- 
ments were performed using a Quantum Design Physical 
Properties Measurement System equipped with a vibrat- 
ing sample magnetometer option. Polycrystalline sam- 
ples were placed into gelatin capsules and affixed to the 
end of a carbon fibre rod. Data were taken on warm- 
ing from 2 K in a fixed magnetic field of 0.1 T. The 
datasets for the two compounds are quite distinct (see 
Fig. S3), but both exhibit the broad hump in Xmoi(T) 
at T J expected for highly anisotropic antiferromag- 



nets and caused by the onset of short-range correlations. 
Fits from 300 K down to temperatures just below the 
hump were made to theoretical models based on the 
low-dimensional Heisenberg Hamiltonian of Equation 1. 
The layered compound is found to be in excellent agree- 
ment with the S = 1/2 Q2D model of Woodward et 
al. [S9] with g = 2.21 and J = 8.10 ± 0.01 K, while the 
best fit for the chain-like material is for the quasi-one- 
dimensional model of Johnston et al. (S10) with g = 2.11 




10 100 
temperature (K) 

FIG. S3. Results of temperature-dependent susceptibility 
measurements. Data points are the powder magnetic sus- 
ceptibilities for both materials, red lines are fits to low- 
dimensional models described in the text. 



and J = 8.58 ± 0.01 K. These values are in good agree- 
ment with those derived from the pulsed field magne- 
tization. We obtained similar results by also fitting to 
the earlier models of Bonner & Fisher |S8| (Q1D) and 
Lines [32] (Q2D). 

Electron-spin resonance. Measurements of the g- 
factor anisotropy were performed in a superconducting 
magnet at the National High Magnetic Field Labora- 
tory in Los Alamos using a cavity perturbation tech- 
nique |S13| . The TE102 mode of a rectangular res- 
onator at 71.5 GHz, and a Millimetre-wave Vector Net- 



work Analyser manufactured by ABmm were used. The 
cavity was mounted on a goniometer and placed inside 
standard flow 4 He and single-shot 3 He cryostats. A tem- 
perature dependence (down to 0.5 K) of the spectrum 
along the principal axes was performed to ensure the re- 
ported angular dependance of the resonance is not influ- 
enced by the long range order, but characteristic of the 
paramagnetic state. For both compounds the data shown 
satisfy this condition. 

The magnetic properties of the two materials studied 
are summarized in Table IS3I 
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